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ABSTRACT 

We investigate roles of Alfvenic waves in the weakly-ionized atmosphere of hot Jupiters by carry¬ 
ing out non-ideal magnetohydrodynamic (MHD) simulations with Ohmic diffusion in one-dimensional 
magnetic flux tubes. Turbulence at the surface excites Alfven waves and they propagate upward to 
drive hot (~ 10^ K) outflows. The magnetic diffusion plays an important role in the dissipation of the 
Alfvenic waves in the weakly ionized atmosphere of hot Jupiters. The mass-loss rate of the sponta¬ 
neously driven planetary wind is considerably reduced, in comparison with that obtained from ideal 
MHD simulations because the Alfvenic waves are severely damped at low altitudes in the atmosphere, 
whereas the wave heating is still important in the heating of the upper atmosphere. Dependence on 
the surface temperature, planetary radius, and velocity dispersion at the surface is also investigated. 
We find an inversion phenomenon of the transmitted wave energy flux; the energy flux carried by 
Alfven waves in the upper atmosphere has a nonmonotonic correlation with the input energy flux 
from the surface in a certain range of the surface temperature because the resistivity is determined 
by the global physical properties of the atmosphere in a complicated manner. We also point out that 
the heating and mass loss are expected only in limited zones if the open magnetic field is confined in 
the limited regions. 

Subject headings: planets and satellites: atmospheres — magnetohydrodynamics (MHD) — planets 
and satellites: magnetic fields — planets and satellites: gaseous planets 


1. INTRODUCTION 

Recently, a large number of exoplanets have been found 
by various detection methods, and it have been learned 
that the exoplanets are rich in diversity. Some of them 
have very small semi-major axes and masses that are sim¬ 
ilar to the Jupiter, and they are known as hot Jupiters. 

Transit observations are a very useful method to know 
not only the radius of planets and their orbital period, 
but also information about their atmosphere. For ex¬ 
ample, transit observations in the UV band provide in¬ 
teresting physical properties of the upper atmosphere. 
Absorption in the H Lya have bee n confirmed in some 
hot J upiters, such a s HD 209458b (fVidal-Madiar et all 
1200, Jfl . HD 18 9733b (iLecavel ier des Etangs et al.ll2010tl . 


and 55 Cnc b (jEhrenreich et al]l2012h . The transit depth 
corresponds to the apparent radius of planets. While 
the transit depth in optical wavelength at HD 209458b 
is 1.5%, that corresponds to ~ 1.38 i?j, the transit 
depth in the H Lyo: line is ^ 15% , that corresponds to 
' 4.3 Rj (|Vidal-Madiar et aI]l2003D. Similar phenomena 


are re ported at HD 189733b ( Lecaveli er des Etanes et"^ 
[2oToh . and a hot Neptune, G.I 436b (|Kulow et al.ll20l5b 


These observations suggest that the existence of the 
extended exosphere that is filled up with the high- 
temperature atomic hydrogen gas, and mass loss takes 
place. Heavier elements such as C, O, Si, and Mg are 
also detect ed in the UV band in the escaping atmo- 
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Although the detailed mechanism of the strong mass 
loss from hot Jupiters is still unclear, various theo- 
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retical models have been developed. A leading ex¬ 
ample of them is that energy-limited escape by X- 
ray and extreme ultraviolet (X UV) from a centra l 
star, so-called XUV-driven escape (|Lammer et al.l[200^ . 
The XUV-driven escape have been studied by v arious 
ways : for e xample, including chemical reactions (jYelle I 
|2004 l2006f) . considering radiation pressure on the es¬ 
caping atmosphere and the charge exchange between 
neutral particles in the escapi ng atmosphere and pro¬ 
tons in the hot stellar wind (IH olmstrbm et al. l l2008l 
Ekeiiback et al.ll20I0HBourrier fc Lecavelier des Etangs I 

20I,3HTremblin fc Chianer Il201,3ll . Additionally, it is suer- 
gested that XUV-driven escape ca n occur not only in hot 
Jupiters but also in cooler planets. IChadney et al.l (|2015f) 
demonstrated that hydrodynamic escape driyen by XUV 
can take place at the gaseous planets that orbit beyond 
I AU in some cases if the central stars are actiye in the 
X-ray. 

Furthermore, physical properties and phenomena re¬ 
lated to magnetic fields of exoplanets haye been inyes- 
tigated in more detail by both theoretical and obserya- 
tional studies. For the theoretic a l work s, for example, 
lAdams I (|201Ill ; I Owen fc Adams I 1)20141 1 suggested out¬ 
flow from the upper atmosphere will be magnetically- 
controlled, thus the outflow becomes asymmetric and a 
mass-loss rate can be smaller about an order of magni¬ 
tude than the yalue without the magnetic fields, that 
is based on the assumption of the spherical symmet¬ 
ric wind. Also, models of the upper atmosphere with 
the planetary magnetic fields are proposed and the ef¬ 
fects on transit depth and the loss of the angular mo- 
mentum by the planetary m agnetic helds were discussed 
(|Trammell et al.ll201lL l2014ll . Recently, magnetic actiyi- 
ties of exoplanets haye been obseryed, and seyeral works 
suggested that the strength of the magnetic fields of hot 
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Jupiters seems to be s maller than that of Jupiter (e.g., 
iKislvakova et alJl2014l) 

We have also investigated magnetically-driven mass 
loss from gaseous planets by extending MHD simulations 
for wi nds from stars with surface convection, e.g. solar 
wind (j Tanaka et al.l[2?)l^ . Alfven waves, which are ex¬ 
cited by turbulence at the surface, transport the energy 
upward. They heat up the upper atomosphere and drive 
outflows. However, in this paper we assumed the ideal 
MHD approximation, which should be modified to take 
into account effects of the weak ionization in the planet 
atmopshere. In this paper, we study non-ideal MHD 
effects by considering the magnetic diffusion in partially- 
ionized atmospheres. 

We explain our concept of the simulations and numer¬ 
ical method, and several formulae for analysis in Section 
[21 Then we show results on the structure of atmopsheres 
and outflows by inspecting roles of the magnetic resistiv¬ 
ity in the propagation and dissipation of Alfven waves 
in Section |3l In Section [H we discuss effects of Ohmic 
dissipation in the atmosphere and effects of geometry of 
the magnetic flux tubes. In Section |5l we summarize the 
results and refer to our future works. 


2. MODEL DESCRIPTION 
2 . 1 . Numerical Method 

In this work we use a numerical simulation code that 
is originally developed for the calculation of th e acce l- 
eration of solar wind ([Suzuki &: Inutsul^ l2005l 120061) . 
At the surf ace of stars, various types of magne tic waves 
are excited (iMatsumoto fc Suzuki Il2012ll2014l) . The ex¬ 
cited waves propagate upward from the surface, then 
heat up the coronal region and drive the s tellar wind 
by va r ious dissipation processes of the waves (iGoldstein I 
I978t iHevvaerts fc PriesiJ Il983t fTerasawa et al.l Il986t 
Kudoh fc Shibata [19991: iM atthaeus et al.| |1999ih It is 
thought that the Alfven wave is especially important for 
the transfer of the energy. This mechanism of the en¬ 
ergy transport from the surface to the upper atmosphere 
can be applied for objects that have their own magnetic 
field. Although details are uncertain, hot Jupiters are 
expected to have magnetic fields, so this mechanism is 
applicable to hot Jupiters. In our previous research, we 
have extended the simulation code to calculate the atmo¬ 
spheric e scape from gaseous p lanets, especially from hot 
Jupiters (I Tanaka et al .11^141) . In addition, this mecha¬ 
nism is applied for the calculation of the structure of the 
atmosphere of brown dwarfs, and it is implied that the 
heating by the dissipation of MHD wav es is also impor¬ 
tant for the brown dwarfs’ atmosphere ([Sorahana et al.l 

IMl). 

As for our MHD simulations we basically follow 
([Tanaka et al.ll201^ . except that we consider magnetic 
diffusion; we solve time-dependent equations and treat 
the propagation and dissipation of the MHD waves in 
a super-radially open magnetic flux tube in the atmo¬ 
sphere. We adopt the same function al form for the super - 
radial expansion factor, /(r), as in (I Tanaka et al.ll2 014l). 
which was originally introduced bv ([Kodd fc Holzer I 
[l97l). 


f{r) 
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where hi is the typical height of closed magnetic loops 
of, and the subscript 0 means the surface. 

The effects of resistivity explicitly appear in the energy 
equation and induction equation: 
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where p, v, p, e, and B are the density, velocity, pres¬ 
sure, specific energy, and magnetic field strength, re¬ 
spectively. The subscripts r and T denote the ra¬ 
dial and perpendicular components; d/dt and d/dt 
denote the Lagrangian and Eulerian derivatives. G 
and M* denote the gravitational constant and the 
mass of a central object, and rj is the magnetic re¬ 
sistivity. Fc and are the thermal conductive 
flux, a nd radiative cooling (iLandini fc Monsignori-Fossi" 


199(11 : [Sutherland fc Dop ib^^ l 199,11 lAnderson fc Athav 
T9M see [Tanaka et al.l ([2014f ) for detail). In these equa¬ 
tions, the effect of the super-radial expansion of the mag¬ 
netic flux tubes appears as the factor r-^. Note that the 
term with p in the energy equation ([2]) indicates Joule 
(Ohmic) heating and the term with r] in the induction 
equation ([ 3 ]) denotes diffusion of magnetic held lines. 

The most dominant origin of the magnetic resistivity in 
the atmosphere of hot Jupiters is the collision between 
electrons and neutral particles. In such circumstances 
the magnetic resistivity 77 depends on gas temperature 
and an ionization degree as follows; 


-q 


234 \/^W (cm^s-i) 


(4) 


where = ndn is the ionization degree (e.g., 
iBlaes fc: Balbus I I1994D . Tie is number density of elec¬ 
trons, and n is that of neutral and ionized hydrogen. If 
the second term on the right-hand side of Equation ([3]) is 
dominant over the hrst term, that is, the magnetic resis¬ 
tivity is large, the magnetic waves will diffuse out. Since 
the gas temperature and ionization degree will be low in 
the atmosphere of planets, the magnetic resistivity might 
not be negligible. Therefore, the resistive MHD calcula¬ 
tion is needed to evaluate the non-ideal effects on the 
structure of the atmosphere and the mass loss. 

Ionization degrees are determine d by a method origi¬ 
nally developed for Betelgeuse by [Hartmann fc Avrett I 
(see also [Harper et al.l[2009[). We ad opt the solar 


abundance gas ([Anders fc Grevesse1ll989t) ^. and calcu¬ 
late the ionization and recombination of H, C, Na, Mg, 
Al, Si, S, Ca, Fe, and K in the gas phase. Since the 
typical surface temperature of hot Jupiters is ~ 1000 K 
and it is low to ionize hydrogen, heavy elements such as 

^ Note: the “so lar composition” is updated recently by 
lAsplund et al.l 1120091') 
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Na and K are the dominant ionizing sources in the lower 
atmosphere. 

The inner boundary of the simulations is set at 
the pho tosphere in the solar an d stell a r wind calcu¬ 
lations (I Suzuki fc Inutsukal l2005l 120061 : ISuzuki I l2007t 
LSuzuki et al.l l2013h . In this paper, we set the inner 
boundary at the position where po = 10®dyncm“^(= 
O.lbar). We fix the temperature at the inner boundary, 
and it is treated as a parameter. The density at the inner 
boundary is determined by the given temperature and po. 
The outer boundary radius is set to 360 planetary radii, 
and the number of the grids is 6000. In the most inner 
region, dr corresponds to 0.0001 planetary radii, and dr 
increases gradually as the distance from the surface. 

We set the planet mass to be the Jupiter mass through¬ 
out this paper. We inject velocity perturbations at 
the inner boundary, which excite upward propagating 
Alfven waves, and their value is treated as a parame¬ 
ter. In most cases, amplitude of the injected pertur¬ 
bations is 20 % of the sound speed at the surface. In 
the case of the sun, the velocity dispersion at the sur¬ 
face is caused by surface convection. On the other 
hand, several previous works have suggested that the 
convective-radi ative boundary lies d eep region of the at- 
mospherefe.g.. [furrows et al.l 1200^1 . However, three- 
dimensional calculations of the atmospheric circulation 
of hot Jupiters suggest the existence of sup ersonic equa¬ 
torial flow(e.g., [S^wma^&i^uniot] [2002|). Therefore, 
it is expected that turbulence can be created at the sur¬ 
face by the strong wind. We assume a broadband spec¬ 
trum of the perturbations in proportion to l/i^, where v 
is frequency of perturbations. To solve the MHD equa¬ 
tions, we adopt a second-or der Godunov sche me with the 
Method of Characteristics (|Sano et al.lll999ll . 

2 .2. Useful Formulae 

We here summarize several useful formula to analyze 
and evaluate simulation results. 


2 .2.1. Scale Height 

We treat the surface temperature Tq, the radius of a 
planet Rp, and the velocity dispersion 5v at the surface 
as parameters. We discuss the dependence of the atmo¬ 
spheric structure on these parameters. To evaluate them, 
we introduce the scale height of the atmosphere. The 
equation of hydrostatic equilibrium of the atmosphere is 


1 dp 
p dr 


GM„ 


= 0 . 


(5) 


Mp is the mass of a planet. If we assume isothermal 
condition, the density profile can be expressed as follows; 


p = po exp 


= Po exp 


r Rp Rp 


Ho r 
mg (r — Rp) Rp 


ksTo 


( 6 ) 


where po is the density at the surface, r = Rp and Hq 
is the pressure scale height of the atmosphere that is 
denoted as 

Ho = (7) 

mg 

where is the Boltzmann constant, m is mean molecu¬ 
lar weight, and g is the gravitational acceleration at the 
planetary surface. 


2 .2.2. Poynting Flux 

Behavior of Alfvenic wave in the atmosphere, such as 
propagation, reflection, and dissipation, is a key to un¬ 
derstand the atmospheric structure. In the ideal MHD 
simulations the dissipation of the Alfvenic waves is done 
by the nonlinear mode conversion. The fluctuations of 
the magnetic pressure associated with the Alfvenic waves 
excite compressive waves. These waves are eventually 
steepen to shockl ets, which finally dis sipate and heat the 
surrounding gas (|Tanaka et al.ll20l1l . This process has 
been examined in terms of the heating of the solar corona 
and the acceleration of the solar wind |H ollweg 1 119821 : 
iKudoh k. Shibata II1999I: ISuzuki fc Inutsuka 1120051) . Re¬ 
cently, a radio observation by JAXA’s spacecraft Akat- 
suki actually reveals radial distribution of compressive 
waves in the solar coro na, which supports this scenario 
(|Mivamoto et al.ll2014f) . 

In the resistive MHD simulations, the magnetic diffu¬ 
sion also contributes to the dissipation of MHD waves. 
When To is lower, the dissipation of MHD waves in the 
low altitude where the magnetic resistivity is high affects 
the structure of the atmosphere. Therefore, an energet¬ 
ics argument of propagating MHD waves in the resistive 
MHD condition is important to understand the proper¬ 
ties of the atmospheric structure and mass loss. 

To evaluate the propagation, reflection, and dissipa¬ 
tion of the Alfvenic wave energy in the atmosphere, we 
summarize the Poynting flux carried by Alfven waves 
below. The net Poynting flux arising from magnetic ten¬ 
sion, which represents the energy flux of Alfven waves 
measured from the comoving frame, can be written as 
follows; 


Tp = — 


Jtt 


( 8 ) 


The Poynting flux can be divided into an outward (par¬ 
allel with Br) part and an inward (antiparallel with Br) 
part (fjacaues I Il97^ ICranmer et al.l [200^ ISuzuki et al.l 
l2013l) . In order to evaluate the inward and outward flux, 
we introduce Elsasser variables, 


z±=v±T 


Bi_ 

\f^P 


(9) 


By using the Elsasser variables, the net Poynting flux is 


-pp = \p (4 - ^-) ^A- 


( 10 ) 


In this equation, pzi^vx/^ and pzfvx!^ correspond to 
the outward and inward components of the Poynting flux, 
respectively. Although in our simulations we inject the 
outward component of Alfvenic waves from the surface, 
the waves are eventually reflected inwardly. 


3. RESULTS 

3.1. Role of Ohmic Dissipation 

In resistive MHD simulations, the propagation and dis¬ 
sipation of MHD waves are affected by magnetic diffu- 
sivity, and accordingly the atmospheric structure would 
be modified. We compare the ideal MHD and resistive 
MHD calculations, and investigate how the magnetic dif¬ 
fusion affects the structure of planetary atmospheres and 
outflows. 
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Fig. 1.— Comparison of the atmospheric structure between 
ideal and resistive MHD calculations, (a) Temperature structure, 
(b) density profile, and (c) radial velocity profile. The solid and 
dashed lines correspond to resistive and ideal MHD cases, respec¬ 
tively. The red and blue lines correspond to the surface tempera¬ 
ture, To =1000K and 2000K, respectively. 

Figure [T] compares the time-averaged atmospheric 
structure of resistive MHD cases with that of ideal MHD 
cases. Note that all atmospheric structures in this pa¬ 
per are time-averaged, and all properties have large time 
variability. The gas density in the upper atmosphere de¬ 
creases more gradually in the ideal MHD cases because 
a larger amount of Alfvenic wave energy can reach the 
upper atmosphere and lifts up the gas by the magnetic 
pressure associated with the waves. On the other hand, 
in the resistive MHD cases Alfvenic waves are damped 
in the cool atmosphere at low altitudes because of the 
high resistivity, and then, the gas density sharply drops 
without the support from the magnetic pressure. The 
Alfvenic waves that enter the upper atmosphere con¬ 
tribute to the acceleration of planetary winds. The wind 
speeds are quite slow in the resistive MHD cases, while 


the wind are accelerated to several tens of kilometer per 
second in the ideal MHD cases. As a result, the mass-loss 
rates in the resistive MHD cases are much smaller than 
those in the ideal MHD cases, because the density in the 
upper atmosphere is around two orders of magnitude less 
than that of the ideal MHD cases, and the velocity is also 
much smaller. If a planet is isolated or have a large dis¬ 
tance from the central star, the planetary wind does not 
stream out in the resistive MHD condition except for in¬ 
termittently driven outflows observed in the simulations 
because the average velocity is smaller than the escape 
velocity. However, if the planet is a hot Jupiter having 
small semi-major axis, the gas overflows easily from the 
small Roche lobe and is blown away by the stellar wind 
from the central star. 

In contrast to the large difference in the p and pro¬ 
files, the temperature structures of the resistive and ideal 
MHD cases are not so different each other. Although the 
locations of the temperature inversion are slightly differ¬ 
ent, the upper atmospheres are heated up to ~ 10^ K in 
both ideal and resistive MHD cases. Therefore, the heat¬ 
ing by the dissipation of MHD waves is still important 
in the atmospheres of the resistive MHD cases. 

3.2. Dependence on Surface Temperature 

We show the relation between the surface temperature, 
Tq, of gaseous planets and the atmospheric structure, and 
then we discuss the dependence of the mass-loss rate. In 
this subsection we fix the radius to the Jupiter radius 
and the amplitude of the velocity dispersion to 20% of 
the sound speed at the surface. Therefore, larger Tq cor¬ 
responds to larger MHD wave energy that is deposited 
to the magnetic flux tube. Tq is thought to be mainly 
determined by the irradiation from the central star, more 
specificly, determined by the combination of the effective 
temperature of the central star and the semi-major axis 
of the planet. This implies that, in the framework of 
the wave-driven planetary winds, when we take plane¬ 
tary systems that have a similar central star, properties 
of the mass loss from hot Jupiters strongly depend on 
the distance from a central star. 

3.2.1. Atmospheric Structure 

The dependence of the atmospheric structure on Tq is 
shown in Figure O In these simulations, the radius and 
mass of the planets are set to the radius and mass of 
Jupiter. In all cases, the gas temperatures in the lower 
atmosphere are almost isothermal, and they rise rapidly 
in the upper atmosphere because of the dissipation of 
the Poynting flux carried by Alfvenic waves. The dissi¬ 
pation of the energy in the upper atmosphere makes a 
high temperature (^10000 K) corona-like region around 
the hot Jupiters. In all the cases, a low-temperature re¬ 
gion forms just below the location of the temperature 
inversion because of the adiabatic expansion of the gas. 
As shown in Figure[2](a), the altitude of the temperature 
inversion is higher for higher Tq. This is mainly caused 
by the difference of the density in the upper atmosphere. 
In cases with higher Tq, the density in the upper atmo¬ 
sphere is higher because of the larger scale height; the 
density decreases more gradually with altitude (Equa¬ 
tion ®). Since larger heating is required to heat up 
the denser gas, the temperature inversion is located at a 
higher altitude in cases with larger Tq. 
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Fig. 2.— To dependence of the atmospheric structure, (a) 
Temperature structure, (b) density profile, and (c) radial velocity 
profile. The horizontal axis denotes the distance from the planetary 
surface that is normalized by Rp in logarithmic scale. The solid, 
dotted, dashed and dot-dashed lines correspond to Tq = 1000 K, 
1600 K, 2000 K, and 2500 K, respectively. 

While the temperature and density panels in Figure [5] 
show more or less similar profiles each other, there is a 
large difference in the radial velocity profiles. The cases 
with low surface temperature Tq < 2000 K yield quite 
slow wind velocities < 3 km s“^. In contrast, in the case 
with high Tq of 2500 K, gas is largely accelerated to at¬ 
tain much faster velocity, ~ 10 km s“^, than the other 
three cases. These results of the velocity and density 
profiles suggest that the wind structure can be catego¬ 
rized into two regimes; while slow and weak planetary 
winds stream out from cooler planets, fast and strong 
planetary winds are accelerated from hotter planets. We 
discuss this difference in Section 13.2.21 

Next, we examine the ionization degree that deter¬ 
mines the magnetic resistivity in the atmosphere. To 
evaluate the non-ideal effects in the weakly-ionized at¬ 


Fig. 3.— To dependence of the (a) ionization degree and (b) 
magnetic Reynolds number in the atmosphere. The horizontal axis 
denotes the distance from the planetary surface that is normalized 
by the Rp in logarithmic scale. The solid, dotted, dashed and dot- 
dashed lines correspond to Tq = 1000 K, 1600 K, 2000 K, and 2500 
K, respectively. 

mosphere, we introduce a magnetic Reynold number, 

= — ( 11 ) 

where v and L are a typical velocity and length in the sys¬ 
tem. We use the scale height of the atmosphere (Equa¬ 
tion dZD) and the amplitude of the velocity dispersion at 
the surface for L and u, respectively. Figure [3] shows the 
dependence of Xe and Rm on Tq. Since the typical value 
of Tq is too low to ionize hydrogen atoms, the ionization 
degree is very low in the lower atmosphere. Xe ~ 10“^*^ 
in the cases with Tq = 1000 K and 1600K, and Xe is still 
10“® even for Tq = 2500 K. In this temperature region, 
the ionization degree is supported by the ionization of 
alkali metals, such as sodium and potassium, that have 
relatively low ionization potential. This low x^ leads to 
high magnetic resistivity and small magnetic Reynolds 
number. As a result, the magnetic Reynolds number in 
the lower atmosphere is very small, Rm in the 

cases with Tq = 1000 K and 1600K, and still Rm 1 
for To = 2500 K, although it is rjIOO times larger owing 
to the higher Xe- The resultant low Rm means that the 
non-ideal effects are important in the lower atmosphere. 
In the upper atmosphere, however, high gas tempera¬ 
ture and low gas density cause a sudden surge of Xe- 
Xe increases gradually in the region below the tempera¬ 
ture inversion region as p decreases with altitude. Xe is 
very high above there because the temperature rises high 
enough to ionize the hydrogen atoms. 
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3.2.2. Two Regimes of the Poynting Flux 

Figure 0] shows the net Poynting flux in the atmosphere 
for cases with different Tq. The net outgoing Poynting 
flux decreases rapidly with an increase of the altitude in 
the lower atmosphere, and it becomes almost constant in 
the upper atmosphere. The Poynting flux injected at the 
surface increases monotonically with Tq, but its behavior 
in the atmosphere is not simple; the net Poynting flux 
in the upper atmosphere is smaller for higher Tq when 
To < 1600K (a), but it is larger for higher Tq when 
To > 1600K (b). The Alfvenic waves are damped in the 
lower atmosphere where the magnetic resistivity is high. 
The scale height of the atmosphere plays a key role in 
the resistive dissipation of the Alfvenic waves. 



Distance from surface ((r-Rp)/Rp) 



Fig. 4.— Net outgoing Poynting flux as a function of the distance 
from the surface for different Tq cases, (a) presents cases with lower 
To < 1600 K and (b) presents cases with higher To > 1600 K. The 
Poynting ffux is multiplied by /rg to separate the effect of 

the adiabatic expansion of the flux tube. 

The scale height is larger for higher Tq , which gives the 
large density in the upper atmosphere. On the one hand 
the larger density in the upper atmosphere leads to a 
larger amount of the mass loss, but on the other hand it 
gives a negative effect on the wind through the enhance¬ 
ment of the resistive dissipation of the Alfvenic waves. 
The ionization degree is lower in the higher-density gas 
because of the efficient recombination, and consequently 
the magnetic resistivity is relatively higher. In Figure 
[3] (b), the magnetic Reynolds number of the case with 
To = 1000 K is low in the lower atmosphere, but it in¬ 
creases quickly in the upper atmosphere. In the case 
with To = 1600 K, it is as low as the Tq =1000 K case in 
the lower atmosphere, but it increases more slowly than 


in the 1000 K case. Therefore, the magnetic Reynolds 
number remains relatively lower in the almost entire re¬ 
gion than in the Tq =1000 K case. Because of the larger 
resistivity, a smaller fraction of the injected wave energy 
reaches the upper atmosphere in the Tq = 1600 K case, 
compared to the Tq = 1000 K case. This is the main 
reason why the transmitted Poynting flux shows the neg¬ 
ative tendency on the surface temperature for Tq < 1600 
K. However, for Tq > 1600 K, the Poynting flux exhibits 
a positive correlation with Tq because the magnetic re¬ 
sistivity is lower for higher Tq. 



Fig. 5.— Comparison of the inward and outward Poynting fiux 
on distance from the surface for the cases with Tq = 1000 K (red), 
1600 K (green), and 2000 K (blue). The Poynting flux is multi¬ 
plied by r^/(r) /tq in order to remove the effect of the adiabatic 
expansion of the flux tube. The solid lines are the outward energy 
flux F+, and the dashed lines are the inward energy flux F —. 

Figure [5] compares the outward and inward energy flux 
the Alfvenic waves leauation ITUl) . Here we show three 
cases with Tq of 1000 K, 1600 K, and 2000 K. Every case 
shows that the inward energy flux tracks the outward 
energy flux with a slightly smaller level in the almost en¬ 
tire region of the atmosphere. The net outgoing Poynting 
flux in Figure 0] is considerably smaller than the outward 
energy flux in Figure [3 These results indicate that the 
injected outgoing Alfven waves are reflected back down¬ 
ward through the propagation in the atmosphere. A tiny 
fraction of the injected Alfven waves can reach the upper 
atmosphere, which contributes to the heating and the ac¬ 
celeration of the gas. In general, a small fraction of the 
injected Poynting flux by Alfven waves is transmitted to 
the upper atmosphere because they suffer dissipation and 
reflection, whereas the resistivity is one of the keys that 
control the dissipation of the waves. The fraction that 
enters the upper atmosphere determines the properties 
of the wind. For instance, if a large amount of the MHD 
wave energy reaches the upper atmosphere, it will rise 
the density there and drive outflows. 

The effect of the resistive dissipation of the 
Alfvenic waves is also seen in Figure [SJ Both dissipation 
and reflection of the waves decrease the Poynting flux. 
The Poynting flux in the case Tq =1000 K (red lines) 
decreases much faster than that of the case Tq =2000 K 
(blue lines), which is caused by both resistive dissipation 
and reflection. The higher magnetic resistivity in the 
1000 K case makes the Alfvenic waves dissipated faster, 
and the steeper density gradient results in strong reflec¬ 
tion. In the 1600 K case (green lines), the Poynting flux 
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quickly drops, which is much faster than in the other two 
cases. Since the scale height is larger than that of the 
To =1000 K case, the reflection of the Alfv&ic waves is 
suppressed. In this case the resistive dissipation of the 
Alfvenic waves is the primary reason of the rapid drop 
of the Poynting flux, and the amount of the energy that 
can reach the upper atmosphere is quite small. 

3.2.3. Time Variability of Mass-loss Rate 

Our simulations show large time variabilities of the at¬ 
mospheric structure. Figure [6] shows the mass-loss rate 
for different Tq cases with time. One can see that in the 
cases with lower Tq < 1600 K the outflow is not time- 
steady, but the atmospheric material infalls during cer¬ 
tain periods, which is seen as discontinuous lines. For ex¬ 
ample, in the case with Tq = 1000 K, the outflow phases 
are relatively shorter and it is comparable to the total 
duration of the infall phases. In the case with Tq = 1600 
K, the integrated duration of the outflow phases is longer, 
and the infall phase disappears for higher Tq > 2000 K. 



Fig. 6.— Time evolution of the mass-loss rates for the cases 
with To = 1000 K (red), 1600 K (green), 2000 K (blue), and 2500 
K (violet). Thick solid lines show the outflow rates, and dashed 
lines show the infall rates. On the left the mass-loss rate is shown 
in units of solar mass per year and on the right in units of gram 
per second. 

These results suggest that in the cases with low Tq < 
1600 K the atmosphere is almost hydrostatic because the 
damping of Alfvenic waves is too severe to drive steady 
wind, and the outflows occur only intermittently. For 
higher Tq the intermittency of the planetary winds dis¬ 
appears, more steady outflows are obtained. 

In Figure [7] we compare the Poynting flux by the 
Alfven waves in the active phase during 23.6-31.6 hour 
and that in the infall phase during 17.8-20.2 hour in Fig- 
ure[6l In the outflow phase, a larger amount of the Poynt¬ 
ing flux is carried by the Alfven waves into the upper 
atmosphere, while in the infall phase, the Poynting flux 
in the upper atmosphere is strongly reduced around four 
orders of magnitude. This result implies that the trans¬ 
mission of the Poynting flux into the upper layer controls 
the on-off nature of the atmospheric escape for the low 
Tq cases. However, we would like to note that the case 
with Tq = 1000 K (not shown) gives more complicated 
behavior; Poynting flux during outflow phases does not 
always exceed Poynting flux during infall phases bucause 



Fig. 7.— Comparison of the Poynting flux in the outflow (solid) 
and infall (dashed) phases for the case with Tq = 1600K. The 
Poynting flux is multiplied by r^f(r) jr^ in order to remove the 
effect of the adiabatic expansion of the flux tube. 

of the rapid temporal variations. 

3.3. Dependence on Planet Radius 

Here we describe the dependence of the atmospheric 
structure and the mass-loss rate on the radius of gaseous 
planets Rp. A number of theoretical works on the ther¬ 
mal evolution of the gaseous planets suggest that radii 
of gaseous planets, which tend to be around Jupiter ra- 
dius, are affected b y surrounding circumstances (e.g., 
iBurrows et al.l 120071) . In addition, observations have 
shown that some hot Jupit ers have enormous ly larger 
radii than that are expected (|Baraffe et al.l[2r)T?il) . There¬ 
fore, we treat Rp as a parameter and investigate the de¬ 
pendence on Rp. The mass-loss rate is very small for low 
Tq, therefore we set Tq =2000 K in this subsection. We 
also fix the amplitude of the velocity dispersion at the 
surface, Sv = 0.2cs. 

Figure |8] shows the dependence of the atmospheric 
structures on Rp. At low altitudes the temperature is 
almost constant (« surface temperature) in all the cases. 
It decreases initially but rises eventually by the heating 
owing to the wave dissipation. The temperature inver¬ 
sion is located at a higher altitude for larger Rp, because 
the scale height is larger; the slower decrease of the den¬ 
sity leads to higher density in the upper region (Figure 
B: which is not heated up to high temperatures. Accord¬ 
ingly, the bending location of the density gradient, which 
corresponds to the location of the temperature inversion, 
is systematically higher for larger Rp by the same reason. 

In all the cases, the temperature of the upper atmo¬ 
sphere is elevated to 10^ K, by the dissipation of MHD 
waves. However, the wind velocity is quite slow. Partic¬ 
ularly in the cases with the large Rp > 1.2Rj, the gas 
in the atmosphere is almost static or even infalls partly, 
whereas the cases with the smaller radius show outflows 
with a few km s“^. 

Figure [9] shows the Rp dependence of the ionization 
degree and magnetic Reynolds number. The ionization 
degree, Xe, is quite low at low altitudes because of the 
low temperature there (~ 2000 K) in all the cases. The 
dominant ion sources in the lower atmosphere are the 
alkali metals. In contrast, Xe increases to 10“^ — 10“^ in 
the upper atmosphere in all cases. In the smaller Rp case, 
Xe starts to increase from a lower altitude and the final 
Xe is also higher in the upper atmosphere, because the 
recombination is slower there due to the lower density. 
































































Fig. 8.— Rp dependence of the atmospheric structure, (a) 
Temperature structure, (b) density profile, and (c) radial velocity 
profile. The narrow solid, solid, dashed, dotted, and dot-dashed 
lines correspond to Rp = O.SRj, l.ORj, 1.2Rj, l.GRj, and 2.0Rj, 
respectively. Distance from the surface (horizontal axis) is normal¬ 
ized by Rp each case. 

The mass-loss rate remains very small in these cases. 
The case with Rp = Rj gives the mass-loss rate of 2.6 x 
lO“^®M 0 yr“^, and the case with 0.8 Rj gives 2.3 x 
10“^® Mq yr“^, that correspond to 1.7 x 10^gs“^, and 
1.4 X 10' gs“^, respectively (here, Mp = Mj and Tq = 
2000 K). In the ideal MHD calculations, the mass-loss 
rate in the case with Rp = R j and Tq = 2000 K is 
1.9 X iQio gs~^ (|Tanaka et al.l 120141) . Our results show 
that the large resistivity as a result of the low ionization 
degree in the low atmosphere significantly suppresses the 
mass loss from hot Jupiters. However, we would like to 
note that the wind shows large fluctuations. There is 
a possibility remained that the planetary wind streams 
out in an intermittent manner. Also, the Roche lobe of 
a hot Jupiters is small because of the small separation 




Fig. 9.— Rp dependence of the (a) ionization degree and (b) 
magnetic Reynolds number in the atmosphere. The narrow solid, 
solid, dashed, dotted, and dot-dashed lines correspond to Rp = 
O.SRj, l.ORj, 1.2Rj, 1.6Rj, and 2.0Rj, respectively. Distance 
from the surface (horizontal axis) is normalized by Rp each case. 

from the central star. Then, once a slow gas material 
overflows the Roche lobe, it will be blown away by the 
stellar wind from the central stars. 

3.4. Dependence on Velocity Dispersion 

The velocity amplitude at the surface, Sv, is also an 
important parameter that controls the structure of the 
atmosphere and the wind of a planet. In this subsection, 
we fix Rp = Rj, Mp = Mj, and Tq = 2000 K. In Figure 
m we show the atmospheric structure of the cases with 
different Sv. Cases with different Sv give similar density 
and temperature profiles. On the other hand, larger Sv 
gives faster wind velocity. Because of the similarity of 
the density and the temperature shown in Figure [TUI the 
ionization degree and the magnetic resistivity are also 
similar in all the cases. Therefore, it can be said that 
the difference of Sv at the surface has small effect on the 
structure of the atmosphere except for the radial velocity. 

Though the difference of the density in the upper at¬ 
mosphere is very small, the wind velocity is faster for 
larger Sv, and it causes the increase of the mass-loss rate 
as shown in Figure [Tl] for cases with Tq = 1000 K and 
2000 K. 

4. DISCUSSION 

4.1. Ohmic Heating in Planetary Atmosphere 

We have investigated the effect of the magnetic diffu¬ 
sion on the dissipation of the Alfvenic waves, which sup¬ 
presses the planetary winds, so far. On the other hand, 
the resistive dissipation also heats up the ambient gas 
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Fig. 10. — Dependence of the atmospheric structure on the 
5v at the surface, (a) Temperature structure, (b) density profile, 
and (c) radial velocity profile. The horizontal axis denotes the 
distance from the planetary surface that is normalized by the Rp 
in logarithmic scale. The solid, dashed, dotted, and dot-dashed 
lines correspond to 5v = 0.2cs, 0.5cs, 0.7cs, and 0.9cs, respectively, 
and Cs is the sound speed at the surface. In all cases, Rp = i?j, 

^1, oiimi(?^laul^i?eatmg, which plays a positive role in 
driving outflows. In this subsection, we examine the role 
of the Ohmic heating in a quantitative manner. 

Current density is 

j = ^ (V X B), (12) 


and electric conductivity a is related with the magnetic 
resistivity rj by 


1 

<j 



(13) 


Then, Ohmic heating rate per unit volume, Qohm = 
is derived as 


Qohm = ^ (V X B)^ 


(14) 



Fig. 11.— Relation between Svfcs and the mass-loss rate for 
To = 1000 K (circle) and 2000 K (filled square). The mass-loss 
rate is shown in units of solar mass per year (left axis) and gram 
per second (right axis). 




Fig. 12.— Comparison of ideal MHD heating rate (red solid), 
cooling rate (green dashed), and Ohmic heating rate (blue solid) 
per unit volume for cases with Tq = 1000 K (a) and 2000 K (b). 
We set 5v = 0.2cs, Mp = Mj, and Rp = Rj. 

Figure [12] shows the Ohmic heating rate in compari¬ 
son with the ideal MHD heating (red solid) and cooling 
(green dashed) rates of the gas. The ideal MHD heat¬ 
ing includes shock heating (see Section 12.2. 2|1 in addi¬ 
tion to adiabatic heating. In the lower atmosphere, the 
Ohmic heating dominates over, or be comparable to the 
ideal MHD heating because of the large resistivity. How¬ 
ever, since the Alfvenic waves dissipate very rapidly, the 
Ohmic heating rate drops drastically as the altitude in¬ 
creases. 
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Fig. 13.— Dependence of the Ohmic heating rate per unit 
volume on Tq. The solid, dotted, dashed, and dot-dashed lines 
correspond to To =1000 K, 1600 K, 2000 K, and 2500 K, respec¬ 
tively. 

Resistivity depends sensitively on temperature. Fig¬ 
ure m compares the Ohmic heating rates for cases with 
different Tq. The higher Tq cases have smaller Ohmic 
heating rate in the lower atmosphere, but it reverses in 
the upper atmosphere, because the dissipation of the 
Alfvenic waves is more gentle at low altitudes in these 
cases and a larger amount of the wave energy reaches the 
upper atmosphere. In addition, as shown in the Figured 
the magnetic Reynolds number is relatively lower in the 
upper atmosphere for higher Tq, which further increases 
the resistive heating. 

However, the Ohmic heating rate is quite small in com¬ 
parison to the shock heating except for the region closed 
to the surface (Fig. IT^ . so this effect is negligible in the 
upper atmosphere and gives little contribution to the at¬ 
mospheric escape. Even in the low atmosphere where 
the Ohmic heating is large, it does not play a primary 
role. The comparison between the ideal and non-ideal 
calculations in Figure [T] (a) shows that the temperature 
profiles in the lower atmosphere are almost same. In 
fact, the temperatures in the lower atmosphere in the 
resistive cases are barely higher than the ideal cases and 
the difference is only several tens of kelvin. In summary, 
the Ohmic heating does not play important roles in the 
entire structure of the atmosphere and the atmospheric 
escape. 


4.2. A Possibility for Bipolar Outflows 

In this work and the previous work, we have fixed the 
shape of the open magnetic flux tube that is described 
in Section EH We presume that the overall structure of 
magnetic field lines of gaseous planets is similar to that 
of the sun. The sun has a multipole magnetic field and 
the propertie s of open magnetic f l ux tubes are o bserved 
in detail fe.g.. lTsuneta et ^l2008Hlto et al.ll20ICll) . How¬ 
ever, the actual structures of magnetic field in exoplan¬ 
ets are unknown. The overall structure of the planetary 
magnetic field affects the shape of the open magnetic flux 
tube. The most critically affected parameter in our mod- 
eli ng is the filling fact or /(r) (see Equation (I) and (2) 
of iTanaka et al.l (j20I4lF . If the planetary magnetic field 
has a simple dipole structure, open magnetic flux tubes 
appear in the magnetic polar regions. In this case, heat¬ 
ing in the upper atmosphere and the acceleration of gas 
flow by MHD wave energy occur only in the polar regions. 
Therefore, the actual appearance of flow is expected to be 


bipolar. In the simple dipole case of the magnetic field, 
the degree of super-radial expansion of the open mag¬ 
netic flux tube should be smaller, and thus the mass-loss 
rate is expected to be smaller because the reflection of 
the energy flux is expected to be enhanced. To inves¬ 
tigate this more quantitatively, a three-dimensional cal¬ 
culation is required, which remains difficult at present. 
On the other hand, the observation of the bipolar outflow 
from an exoplanet is, in principle, possible in the detailed 
analysis of the differential spectroscopy for transiting hot 
Jupiters, unless the magnetic dipole axis is perpendicular 
to both the orbital plane of the exoplanet and the line- 
of-sight. This is because a blue-shifted outflow before 
ingress and a red-shifted outflow after the emersion, or 
vice versa, can be seen in this case. Once the mass flux 
of the bipolar outflow from the exoplanet is observed, 
we will possibly constrain the property of magnetic field 
of the exoplanet by assuming the driving mechanism as 
described in this work. 

4.3. Limitations in the Present Model 

As we mentioned in the introduction, high energy ra¬ 
diation such as X-ray and extreme ultraviolet from a 
central star is important for the heating of upper atmo¬ 
sphere and atmospheric escape from hot Jupiters. Here 
we compare the height where the heating by MHD waves 
is important and the penetration depth of stellar XUV. 
It is useful to introduce Rxuv, the radius at which XUV 
is absorbed. i?xuv means that the optical depth for 
XUV reaches unity at that radius, and it corresponds 
to the radius where the colum n density of hydrogen be¬ 
comes ^ 5 >< fO^J cm“^(e.g., iMurrav-Clav et ah 120091 : 
[Rogers et al.ll20llh . We calculate i?xuv in each calcu¬ 
lations and compare it with the region where the heating 
by the dissipation of the MHD waves becomes important. 



Fig. 14.— Dependence of Rxw oii Fq and comparison with 
Rlow Rsoook- The black squares with the solid line corre¬ 
sponds to /?xuv? the circles with the dashed line and the triangles 
with the dotted lines correspond to Riow Rsoook? respectively. 
We set 5v = 0.2cs, Rp = Rj, and Mp = Mj for this comparison. 


Figure [Til shows Tq dependence of i?xuv and compar¬ 
ison with other considerable radii. As shown in Figure 
Ha), the gas temperature slightly drops due to adiabatic 
expansion and then heated drastically to ^ 10"*^ K. Riow 
is the radius at which the temperature become lowest 
and heating by MHD waves starts to become important. 
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and -Rsoook is the radius at which the temperature be¬ 
come higher than 8000 K. In this range of To, i?xuv 
always locates only slightly above Riow and the value 
varies between 1.05 ^ 1.20i?p. This is consistent with 
previous works on the XU V heating in the atmosp here 
of hot gaseous planets ie.g.. lMnrrav-Clav et al.l[20(M . In 
addition, TgoooK always lies much higher altitude than 
i?xuv- 

These result suggest that the region heated to 10^ 
K by the dissipation of the MHD waves is optically thin 
to XUV photons, and XUV can reach the altitude that 
the temperature becomes lowest in our model. Photo¬ 
ionization, cooling from metal lines, Lyman a play im¬ 
portant roles in the upper atmosphere of hot Jupiters. At 
this stage, our model does not include the effects of stellar 
XUV irradiation, the inclusion of which should improve 
our model. In the present modeling, the results cannot 
explain the observational results of the mass-loss rates 
adequately, since the mass-loss rates are reduced largely 
when we take into account the effects of the magnetic re¬ 
sistivity. This seems to indicate that the resulting mass- 
loss rate strongly depends on the details of the input 
physics. More realistic calculation is critically needed, 
and will be out future work. 

5. CONCLUSION 

We have studied the structure of the atmosphere and 
the wind of hot Jupiters by non-ideal MHD simulations, 
particularly focusing on the effects of the magnetic diffu¬ 
sion. The Alfvenic waves generated from the surface are 
strongly damped via resistive dissipation in the weakly- 
ionized atmosphere at low altitudes and only a tiny frac¬ 
tion of the initial energy can reach the upper atmosphere. 
Consequently, the mass-loss rate is reduced significantly, 
compared to that obtained from the ideal MHD simula¬ 
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tions. However, the Alfvenic waves that survive to trans¬ 
mit into the upper atmosphere still support the temper¬ 
ature inversion and heat up to ~ 10^ K by the wave 
dissipation; the Alfvenic waves are still important in de¬ 
termining the temperature profile even though the mag¬ 
netic diffusion is taken into account. We would like to 
note that the wave heating here is not owing to Ohmic 
heating but the consequence of the shock dissipation of 
the compressive waves that are nonlinearly excited from 
the Alfvenic waves. 

We also investigated the dependence of atmospheric 
properties on the surface temperature, the planet radius, 
and the velocity dispersion at the surface. We found a 
nonmonotonic dependence of the Poynting flux on the 
surface temperature (Figure [U, because the resistivity, 
which controls the dissipation of the Alfvenic waves, is 
determined globally the density and temperature of the 
atmosphere. We point out that the heating and acceler¬ 
ation of gas will take place in magnetic polar regions if 
the planetary magnetic field is simple dipole structure, 
and non-spherically symmetric bipolar outflow will takes 
place. 

In our calculations, the treatment of the radiative cool¬ 
ing is simplified es pecially for low temp erature gas (see 
discussion part of iTanaka et al.l (|2014IH . Therefore, a 
more accurate treatment is needed for future investiga¬ 
tion and now we are going to improve this part. 
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